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Introduction 


An  increased  interest  over  the  past  few  years  in  the  use  of 
ultrafiltration  processes  for  the  treatment  of  surface  and  ground 
water  has  prompted  a  number  of  studies  to  evaluate  economic 
advantages  of  UF  over  other  types  of  treatment  processes.  The 
primary  obstacle  in  making  an  accurate  economic  evaluation  has  been 
the  ability  to  accurately  model  the  response  of  the  UF  process  to 
membrane  fouling.  We  will  interest  ourselves  specifically  with 
hollow  fiber  membrane  modules,  and  attempt  to  accurately  model  UF 
system  parameters  with  a  hydraulic  equivalent  of  the  UF  system. 

Models  based  on  mass  transfer  (film)  theory  and  the  series 
resistance  model  have  been  somewhat  effective,  but  have  had 
specific  limitations  in  their  application  (Cheryan  1986,  Clark 
1991,  Laine  1989).  The  mass  transfer  model  presented  by  Cheryan 
(pp. 84-89)  assumes  system  operation  in  a  region  where  permeate  flux 
is  independent  of  system  pressure.  In  this  model,  the  permeate 
flux  is  described  in  terms  of  boundary  layer  concentration,  gel 
layer  concentration,  and  mass  transfer  coefficient;  although 
possible  to  derive  a  solution  based  on  these  parameters,  the 
solution  technique  is  somewhat  cumbersome  and  does  not  relate 
permeate  flow  to  other  system  parameters,  nor  does  it  allow  for  a 
solution  when  the  system  is  operating  within  the  pressure 
controlled  region.  Because  the  hollow  fiber  membrane  typically 
operates  at  such  pressures,  we  will  not  pursue  further 
investigation  of  the  mass  transfer  model.  The  series  resistance 
model  discussed  by  Clark  and  Laine  has  the  advantage  of  effectively 
modeling  system  response  in  the  pressure  controlled  region,  but 
does  not  allow  for  the  possible  interdependence  of  the  various 
resistance  terms  (Laine,  1989,  p. 66),  nor  does  it  allow  for  non- 
linear relationships  between  flow  and  pressure  within  the  system 
(Clark,  1991,  p. 14) . 

The  motivation  for  this  report  was  a  suggestion  by  Clark 
(1991)  that  simple  membrane  fouling  models  could  be  incorporated  in 
an  overall  UF  system  hydraulic  model.  The  hydraulic  model  we  will 
develop  will  differ  from  the  models  described  above  in  that  we  will 
attempt  to  describe  the  UF  system  as  a  whole  as  opposed  to  modeling 
only  the  membrane  permeate  flux.  As  we  shall  see,  this  approach 
will  allow  for  a  relatively  simple  economic  analysis  of  UF  pilot 
system  operations.  The  hydraulic  model  will  also  offer  the 
advantages  of  allowing  for  the  interdependence  of  the  various  head 
loss  terms  and  for  the  non-linearity  of  head  loss  as  a  function  of 
operating  pressure  and  flow  rate.  The  model  further  allows  the 
user  to  input  any  number  of  membrane  fouling  schemes  directly  into 
the  model  by  varying  a  single  parameter. 

The  hydraulic  model  developed  here  is  based  on  the  application 
of  fundamentals  of  fluid  mechanics  to  the  UF  pilot  system.  The 
head  loss  through  the  membrane  fibers  as  well  as  the  head  loss  of 
the  permeate  passing  through  the  membrane  are  developed  such  that 
they  can  be  expressed  in  terms  of  easily  measured  system  operating 
parameters.    Head  loss  terms  for  system  valves  and  pipes  are 


developed  in  the  same  manner.  A  physical  property  of  the  membrane, 
referred  to  as  the  membrane  constant,  is  introduced,  referring  to 
the  inherent  resistance  of  the  membrane  to  permeate  flow.  Although 
a  physical  constant  of  the  membrane,  it  may  be  regarded  as  a 
property  of  the  membrane  which  one  can  vary  with  either  time  or 
permeate  volume,  thus  allowing  for  the  modeling  of  system  response 
to  various  membrane  fouling  scenarios. 

Once  the  model  has  been  solved,  the  data  output  may  be  used  to 
determine  an  optimal  operating  period  for  the  ultrafiltration  and 
back  flush  cycles.  This  economic  optimization  will  determine  the 
cycle  time  for  which  permeate  flow  per  kilowatt  hour  of  energy 
consumed  is  maximized. 


The  Hydraulic  Model 

Figure  (1)  shows  the  laboratory  set  up  of  the  ultrafiltration 
system  to  be  modeled.  Figure  (2)  is  the  hydraulic  equivalent  of 
the  system.  For  analysis  purposes,  we  shall  consider  a  closed 
system  where  both  permeate  and  retentate  are  returned  to  the  pump 
suction.  We  shall  describe  each  component  of  the  system  in  terms 
of  flow  and  head  loss,  write  the  associated  head  loss  equations, 
and  then  develop  a  solution  procedure  for  several  scenarios. 

System  Components; 

Pump  head:  The  head  delivered  by  the  pump  is  described  by  the 
pump  curve.  Although  the  pump  curve  is  different  for  each  pump,  it 
can  generally  be  described  by  the  following  equation: 

APp  =  aqI+bqc  +  c  CD 

where  AP  is  the  pump  discharge  pressure,  Qt  is  the  flow  generated 
by  the  pump,  and  the  constants  A,  B,  and  C  are  the  results  of  a 
"least  squares"  fit  of  the  pump  curve  (for  the  pump  used  in  this 
model,  A=-36.943,  B=61. 66335,  and  C=63. 40662)  (Clark,  1991). 

Height  differential  head  loss:  The  model  as  constructed  in 
the  lab  had  a  short  height  differential  between  the  pump  and  the 
inlet  to  the  membrane.  The  head  loss  associated  with  the  change  in 
height  is  given  by 

APr  =  pgz  (2) 

where  p  is  the  fluid  density  (kg/m3)  ,  g  is  the  acceleration  due  to 
gravity,  hf  is  the  head  loss  due  to  the  height  differential,  and  z 
is  the  height  differential. 

System  tubing  head  loss:  Head  loss  in  the  system  tubing  is 
described  by  the  Darcy-Weisbach  Equation: 

hf-^  (3) 

1        2gd 

where  L  is  the  length  of  the  tube,  f  is  the  friction  factor,  V  is 
the  velocity  of  the  fluid  through  the  tube,  and  d  is  the  tube 
diameter.   Substituting  APt=pghf  into  equation  (3) ,  we  get 

APt  =  pfLV2  (4) 

c     2d 


recyc le  loop 


permeate 


membrane 


retentote 


pump 


X  backpressure  valve 


recyc le  tank 


Source:    Clark,    1991 
Figure    (1) 


Qt 


AP. 


A  P. 


Pump 


Qr 


A  P. 


AP. 


AP 


bP 


AP, 


Figure    (2 


The  flow  through  the  system  tubing  is  turbulent  (Re~4000)  and 
the  friction  factor,  f,  can  be  determined  through  use  of  the  Moody 
diagram.  For  the  analytical  approach  we  will  use  in  solving  the 
hydraulic  model,  the  Moody  diagram  can  be  approximated  by  the 
following  equation  (White,  p.  313-314) : 


H-'-HM-teaiT 


where  e  is  the  tube  roughness.  We  shall  assume  a  fairly  "smooth" 
tube,  with  a  roughness  of  .0015  mm.  We  shall  further  assume,  for 
purposes  of  this  analysis,  that  the  length  of  the  tube  from  the 
permeate  exit  to  the  pump  suction  is  zero;  this  assumption  reflects 
the  laboratory  condition  of  the  permeate  discharging  directly  into 
a  separate  container. 

Valve  head  loss;  There  are  two  valves  in  the  piping  network; 
one  which  provides  back  pressure  to  the  fiber  module,  and  one  which 
provides  for  a  bypass  of  the  fiber  module.  The  latter  is  included 
in  the  system  to  allow  for  the  use  of  pumps  of  varying  size  in  the 
system.  If  the  pump  is  properly  sized  for  the  membrane  in  use,  the 
bypass  line  is  not  required.  The  head  loss  for  each  valve  is 
described  by  the  following  equation: 

f        2g 

Substituting  AP=pghf   into  equation    (6) ,    we  get 

AP=M¥l  (7) 

2 

In  each  equation,  K  represents  a  proportionality  constant  which, 
unfortunately,  cannot  be  expressed  in  terms  of  roughness  and 
Reynold's  Number  (White,  p. 333).  Values  of  K  are  listed  in  most 
Fluid  Mechanics  references  for  valves  of  various  designs  and  sizes. 
The  valves  we  used  were  not  listed  in  these  references,  and  had  to 
be  measured  experimentally  by  measuring  the  flow  rate  and  the 
associated  pressure  drop  across  the  valve. 

Fiber  Head  loss:  A  'typical'  fiber  is  shown  in  figure  (4). 
The  flow  through  each  fiber  of  the  membrane  is  in  the  transitional 
region  between  laminar  and  turbulent  flow.  The  head  loss  is  again 
described  by  the  Darcy-Weisbach  equation: 

hf  =  lh£  (8) 

f        2gdf 
Substituting   APf=pghf   into  equation    (8) ,    we  get 

2d 


LPf  =    PfL^  (9) 
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For  laminar  flow,  the  friction  factor  is  given  by  the  expression 

5 


f  =  64/R  .  For  turbulent  flow,  the  friction  factor  is  given  by 
equation  (5) .  Again,  we  are  limited  by  not  knowing  the  roughness 
of  the  membrane  fiber.  We  can,  however,  solve  for  roughness  using 
experimental  results.  By  obtaining  measurements  of  membrane  flow, 
we  can  calculate  fiber  velocity  with  the  equation 


*•  ■  i^H-h) 


which  is  simply  the  average  flow  through  one  fiber  divided  by  the 
cross  sectional  area  of  that  fiber;  n  is  the  number  of  fibers  in 
the  membrane  cartridge  (the  membrane  we  used  had  n  =  20  fibers)  ,  Qm 
is  the  inlet  flow  to  the  membrane,  Q  is  the  retentate  flow  rate, 
and  af  is  the  cross  sectional  area  of  a  single  fiber.  Once  fiber 
velocity  is  known,  we  can  use  equation  (9)  to  solve  for  the 
friction  factor: 


fm    2dfAPf 
9Ltv\ 


(11) 


Now  having  both  the  friction  factor  and  the  fiber  velocity,  we 
can  solve  for  e  in  equation  (5)  .  The  results  of  these  calculations 
for  two  different  experiments  are  shown  in  table  (1)  ,  giving  an 
average    e    of   7.64xl0"5  meters. 

TABLE    (1) 
EVALUATION   OF   FRICTION   FACTOR   FROM  VARIOUS   DATA: 

CONSTANTS:  p  df  /i  Lf 

(kg/mA3)         (m)         (kg/m-sec)         (m) 
997.4     0.00093      9.42E-04  1.03 

EXPERIMENT  #   1: 
0*  QP  Qr  Vf  APf  Re  f  e 

(1/min)   (1/min)  (1/min)   (m/sec)   (N/m2) (m) 

1.547    0.387   1.160   0.949   57918.0   934.39   0.117  9.09E-05 

1.621    0.381   1.240   0.994   66881.5   978.97   0.123  1.03E-04 

1.562    0.402   1.160   0.958   59297.0   943.39  0.117  9.21E-05 

1.573    0.393   1.180   0.965   60676.0   950.03  0.118  9.42E-05 

AVG:  9.51E-05 

Source:   Heneghan,  1991 
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The  Moody  Diagram  is  a  plot  of  f  vs.  Re  for  various  e/d  ratios 
(figure  (3) .  Note  that  the  data  plotted  for  the  two  experiments 
lie  in  the  transition  region  between  laminar  and  turbulent  flows, 
though  each  experiment  tends  lie  closer  to  the  laminar  flow  region. 
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Source:   White,  p.  313. 
Figure  (3) 

Permeate  flow  head  loss:  Permeate  flow  through  each  membrane 
pore  is  considered  to  be  laminar  flow  because  of  low  Reynold's 
numbers  (<1800)  obtained  from  experimental  data  for  various 
membranes  (Cheryan,  p. 76-77).  Figure  (4)  shows  a  side  view  of  a 
•typical'  fiber. 


Qm=Qp+Qr 


Figure  (4) 


In  general  terms,  the  head  loss  through  a  single  pore  is  given 
by  the  Darcy-Weisbach  equation: 

hf  =    fAxVl  (12) 

where  f  is  the  friction  factor,  Ax  is  the  membrane  "skin" 
thickness,  V  is  the  velocity  of  the  fluid  through  the  pore,  g  is 
the  acceleration  due  to  gravity,  and  d  is  the  diameter  of  the 
pore.  For  laminar  flow,  the  friction  factor  is  given  by  f  =  64/Re. 
Substituting  into  equation  (12)  gives 

h   sa^  (1J) 

P9dp 

Our  goal  is  to  express  the  permeate  head  loss  in 
experimentally  measurable  terms.  To  do  so,  we  will  express 
equation  (13)  in  terms  of  pore  flow  rate,  pore  area,  and  pore 
radius: 

VP  -  ^  *>-  **P  CIA) 

AP 

Substituting  equation  (14)  into  equation  (13) ,  we  get 

P9aprp 

We  shall  assume  that  the  membrane  in  question  has  a  total  of  "n" 
pores.   The  total  flow  through  the  membrane,  then,  is  given  by 

o,-i«..  <16) 


The  surface  porosity  of  the  membrane  is  defined  as  the  ratio 
of  total  pore  area  to  total  membrane  surface  area.  If  porosity  and 
surface  area  (A)  are  known,  then  the  total  pore  area  may  found: 


A,  .±.„.  *A  (") 


We  now  substitute  equations  (16)  and  (17)  into  equation  (14)  to 
describe  the  total  flow  through  the  membrane: 
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These  expressions  may  now  be  used  to  express  hf  in  terms  of  the 
membrane  area,  porosity,  and  total  permeate  flow: 

(19) 


hf  = 

8nAx£p 

8nAxgp 

pgip€A 

P9apip 

The  permeate  head  loss  across  the  membrane  surface  may  be 
expressed  in  terms  of  a  change  in  pressure  across  the  membrane. 
Because  the  pressure  drops  along  the  length  of  the  membrane  fiber 
(due  to  fiber  head  loss)  ,  we  will  use  the  average  pressure  drop 
across  the  membrane,  hereafter  referred  to  as  the  "trans-membrane" 
pressure: 

AP    =  Illfl  (20) 

m  2 

where  P1  is  the  membrane  cartridge  inlet  pressure,  and  P2  is  the 
outlet  pressure.  The  relationship  between  pressure  and  head  loss 
is  given  by 

APa=pghf  (21) 

Substituting  into  equation  (19) ,  we  obtain 

AP.  -  **¥£*  (22) 

ZpCA 

The  limitation  of  equation  (22)  is  that  membrane  thickness, 
porosity,  and  pore  radius  are  generally  not  known.  These  terms, 
however,  do  represent  a  physical  constant  of  the  clean  membrane.  If 
trans-membrane  pressure  and  permeate  flow  can  be  measured 
experimentally,  we  can  reasonably  estimate  the  unknown  parameters. 
We  shall  then  define  the  membrane  constant,  C,  as 

c  =  -^4  (23) 

erp2 

with  units  of  1/L.   Equation  (22)  then  becomes 

AP  .  *»CQ>  (24) 


Table  (2)  shows  the  calculation  of  the  membrane  constant  for  the 
same  membrane  using  data  from  two  different  experiments  (Heneghan, 
1991,  Adham,  1991)+.  The  results  compare  favorably  with  each 
other,  as  well  as  with  values  obtained  for  various  membranes  from 
manufacturer's  data    (Table    (3)). 


+  The  membrane  used  in  these  experiments  was  BCDA  cellulose  acetate,  with 
a  Molecular  Weight  cutoff  around  100k  daltons,  manufactured  by  Yonaise  des  Eaux, 
Paris,    France. 


TABLE  2 

EVALUATION  OF  MEMBRANE 

CONSTANT 

EXPERIMENT  #  1 

TEST  # 

APt 

°r 

°r 

C 

(N/M2) 

Cl/min) 

CM3/sec) 

1/M 

1 

1.307E+05 

0.3978 

6.630E-06 

1.569E+11 

2 

1.517E+05 

0.3867 

6.445E-06 

1.874E+11 

3 

1.527E+05 

0.3805 

6.342E-06 

1.918E+11 

4 

1.531E+05 

0.4016 

6.693E-06 

1.821E+11 

5 

1.565E+05 

0.3926 

6.543E-06 

1.905E+11 

AVG: 

1.489E+05 

0.39184 

6.531E-06 

1.880E+11 

EXPERIMENT  #  2 

l_ 

TEST  # 

APt 

°r 

°r 

C 

(N/M2) 

Cl/min) 

(M3/sec) 

1/M 

1 

6.647E+04 

0.1450 

2.4167E-06 

2.112E+11 

2 

9.901E+04 

0.2132 

3.5528E-06 

2.141E+11 

3 

8.115E+04 

0.1785 

2.975E-06 

2.095E+11 

4 

5.316E+04 

0.1258 

2.0972E-06 

1.947E+11 

5 

6.757E+04 

0.1528 

2.5472E-06 

2.037E+11 

6 

6.481E+04 

0.1490 

2.4833E-06 

2.005E+11 

7 

1.350E+05 

0.2880 

4.8E-06 

2.160E+11 

8 

3.199E+04 

0.0782 

1.3028E-06 

1.886E+11 

AVG: 

7.490E+04 

0.1663 

2.772E-06 

2.048E+11 

AVERAGE  OF 

THE  TWO: 

1.964E+11 

TABLE  3 


EVALUATION  OF  MEMBRANE 

CONSTANTS 

FROM  MANUFACTURER'S 

DATA: 

PORE 

d 

r 

AREA 

DENSITY 

POROSITY 

Ax 

c 

Cum) 

Cum) 

(um2) 

(PORES /cm2: 

l 

Cum) 

0.9 

0.45 

0.63617251 

20000000 

0.1272345 

11.5 

4 

.463E+08 

0.18 

0.09 

0.0254469 

300000000 

0.0763407 

12 

1 

.941E+10 

0.09 

0.045 

0.00636173 

300000000 

0.01908518 

5.3 

1 

.371E+11 

0.072 

0.036 

0.0040715 

300000000 

0.01221451 

5.4 

3 

.411E+11 

0.045 

0.0225 

0.00159043 

600000000 

0.00954259 

5.4 

1 

.118E+12 

0.027 

0.0135 

0.00057256 

600000000 

0.00343533 

5.4 

Source:   ( 

8 

^heryan , 

.625E+12 

1986,  p. 60 

Solution  of  the  Hydraulic  Model: 

We  now  have  all  elements  of  the  hydraulic  model  described  in 
terms  of  pressure  drop  and  flow  rate,  specifically  the  head  loss 
terms  for  the  pump,  height  differential,  system  tubing,  system 
valves,  the  membrane  fibers,  and  the  membrane  pores  (equations 
1,2,4,7,9,  &  24).  All  that  is  left  to  do  is  solve  for  the  unknown 
pressures  and  flow  rates.  We  shall  examine  four  scenarios:  Laminar 
and  turbulent  flow  solutions  without  use  of  a  bypass  line,  and 
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laminar  and  turbulent  flow  solutions  for  the  case  when  the  bypass 
line  is  used.  Each  scenario  has  been  solved  using  an  iterative 
technique  along  with  a  "QUATTRO  PRO"  spreadsheet  program.  We  shall 
write  the  equations  here,  develop  the  iterative  solution  approach, 
and  then  show  the  results  obtained  through  use  of  the  spreadsheet. 
Note  that  all  pressures  are  measured  in  PSIG,  lengths  in  meters, 
areas  in  meters  squared,  and  all  flow  rates  are  measured  in  liters 
per  minute.   Appropriate  conversion  factors  are  indicated. 

Case  #1;   Laminar  flow,  no  bypass  valve. 

Figure  (2)  shows  the  hydraulic  model  to  be  solved.  The 
unknowns  to  be  solved  for  are  the  membrane  cartridge  inlet  pressure 
P,,  Qt,  Qr,  and  Q  .  We  need,  then,  four  equations  to  solve  for  the 
four  unknowns. 

The  first  unknown  to  solve  for  is  P1 ,  the  inlet  pressure  to 
the  membrane.  P1  may  be  found  by  subtracting  both  the  inlet  tubing 
head  loss  and  the  height  differential  head  loss  form  the  pump  head 
generated  by  the  pump: 

Px  =  AQ2t  +  B£>t  +  C-APz-AP£/       where 

AP*  =  PgZ  (25) 

C  2dca2t 


The  second  unknown,  permeate  flow  (Q  )  may  be  found  by 
substituting  equation  (24)  into  equation  (20) : 

AP  -  Pi  +  P*  -    8>ACgp  (26) 

2       A 


P2  is  simply  the  change  in  pressure  across  the  back  pressure  valve: 

P2  *   AP^  =  ±*!*2L  (27) 

2abp 

Note  that  we  have  substituted  Vr  =  Qr/a^,  where  a^  is  the  cross 
sectional  area  of  the  back  pressure  valve.  This  substitution  will 
make  the  ensuing  calculations  somewhat  simpler.  Retentate  flow, 
Qr,  is  used  because  only  retentate  fluid  flows  through  the  back 
pressure  valve.  We  now  solve  equation  (26)  for  P.,  and  substitute 
equation  (27)  for  P2: 

8\iCQp)         pKbpQ2r 


a- 


2  a2^ 


We  shall  now  introduce  constants  Z  and  Y: 
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p,  =  2  °^-vp\  -    »****'  (28) 


z  =  /ecuy    i    V     i 

\      A     /\  6895  A  60000 


y  = 


fp^V     l    \f      l      ) 

k2A^Pp895A600002J 


(29) 


where  the  term  (1/6895)  is  a  conversion  factor  from  newtons  per 
meter  squared  to  PSI  and  the  term  (1/60000)  is  a  conversion  factor 
from  cubic  meters  per  second  to  liters  per  minute.  Substituting  Z 
and  Y  into  equation  (28) ,  we  get 


P1   =  2ZQ-YQ'Z 


(30) 


The  third  unknown,  retentate  flow,  may  also  be  described  in 


terms  of  P, : 


P1   =  APf  +  APi 


bP 


(31) 


substituting  equations  (9)  and  (27) ,  and  making  the  substitution 
for  the  constant  Y  as  before,  we  get 


64. 
K*. 


2df 


*Y0l 


(32) 


We  now  substitute  R    =   pvd//i: 


Pi  = 


6*H  ) 


pLfVJ 


2d 


f    ) 


*yqI 


(33) 


Substituting  flow  rate  for  velocity  using  equation  (10)  (recall 
n=20)  gives 


Pi  = 


'   -8nL/ 


(Qt*Qr)  +  YQi 


(34) 


Note  that  in  the  case  of  no  bypass  valve,  membrane  flow  Qm  is 
exactly  equal  to  total  flow  Qt.   Defining  a  constant  X: 

(  .8nLfV 

|2 


X  = 


I  d't*t  ) 


\6895/\  60000/ 


(35) 


yields  the  final  form  of  the  equation  for  P1 : 

P1   =  X{Qt*Qz)  +YQJ 


(36) 


Finally,  we  use  the  continuity  equation  to  solve  for  the 
fourth  unknown,  Qt: 
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Ot   =  QP+QT  (37> 

The  four  unknowns  are  now  solved  with  equations  (25, 30, 36, &  37)  in 
the  following  manner:  First,  guess  a  value  for  Qt ,  and  generate  a 
value  of  change  in  pressure  using  equation  (25) .  Then,  solve 
equations  (30)  and  (36)  for  the  quantity  Y*Qr2.  To  do  this,  set 
equations  (30)  and  (36)  equal  to  each  other,  and  substitute  Q  =Q  -Qr 
in  equation  (30) .  We  arrive,  then,  at  the  following  solution  for 

Or* 

_  2P1-QeU  +  2Z)  (38) 

1  X-2Z 

We  can  now  solve  equation  (38)  based  on  the  results  of  the  assumed 
value  of  Qt.  We  then  use  the  continuity  equation  (equation  (37)  to 
solve  for  Q  : 

Qp  =  Qt-Qt  (39) 

We  check  our  estimate  of  Qt  by  solving  for  change  in  pressure  in 
equation  (30) ,  and  iterate  on  Qt  until  the  values  of  P1  given  by 
equations  (25)  and  (30)  match.  These  calculations  are  done  on  the 
QUATTRO  PRO  spreadsheet. 

Case  #  2 :  No  bypass  valve.  Turbulent  flow  through  membrane 
fibers. 

Because  the  flow  through  the  membrane  fibers  is  now  to  be  evaluated 
as  turbulent,  the  equation  for  the  head  loss  through  the  membrane 
fibers  is  given  by 


1     1  2d,  J 


*YQ*  <A°> 


where  the  friction  factor  f  is 


'•hHfclwT 


We  now  define  the  constant  W: 


w  =  -=■ 


i(±£fV_M  («) 


2  I    dt   \\  6  89  5 


Substituting   into  equation    (40) ,    we  get 

p,  =  WfVJ  +  YQ?  (43> 


where 


Vf  -  [°£l°L\       1       \  (44) 

f      {   40a,  j\  60000/ 
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To     solve     the     problem,     we     first     substitute     Q      =Qt-Qr     into 
equation    (30) ,    and   solve    for  Qr.      The   resulting   equation    is 


Or 


-  EMI 


+  4 


2Z£>   -AP 


(45) 


Once  again,  we  guess  a  value  of  Qt  and  solve  for  change  in  pressure 
in  equation  (25) .  Plugging  the  results  into  equation  (45)  will 
give  two  possible  values  of  Qr,  only  one  of  which  is  feasible  (the 
other  will  give  a  negative  value).  We  then  solve  equation  (43)  for 
change  in  pressure  to  check  our  original  guess,  and  then  iterate 
until  the  two  values  match.  Once  again,  the  calculations  are  done 
on  the  spreadsheet  program. 

Case  #  3:   Bypass  valve  in  place,  laminar  flow  solution. 

The  solution  technique  for  this  case  is  similar  to  that  of  case 
number  one;  in  fact,  equations  (25)  &  (30)  still  apply.  Equation 
(36)  remains  the  same,  with  one  exception:  total  flow  leaving  the 
pump  is  now  split  between  the  bypass  valve  and  the  membrane 
cartridge.   The  continuity  equation  now  becomes 

Qt   =  Qe  +  Qm      .where       Qm   =  Qp  +  Qz  (46) 


Equation  (36)  now  becomes 


Pa  =  X{Q^Qr)  +YQ2r 


(47) 


The  addition  of  the  bypass  valve  adds  a  fifth  unknown  to  the 
problem,  the  bypass  flow  (Qs)  .  This  flow  can  also  be  expressed  as 
a  function  of  change  in  pressure  and  the  head  loss  due  to  the 
bypass  tubing  and  the  bypass  valve: 


/        -2\ 


2d, 


PKsVs 


(48) 


Substituting  V  = 


Qs/as, 


Px    =    ', 


!i2dea2B 


2a\ 


(49) 


We  now  introduce  constants  F  &  G: 


-(Hi) 


;00002  J(  6895  J 


(50) 
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G  = 


PL, 


2dBa2s, 


00002j\6895j 


(51) 


Substituting  equations  (50)  and  (51)  into  equation  (49)  yields 

P1   ■  Ql(F+fmG)  <52) 

To  solve  the  equations,  we  again  guess  a  value  of  Qt  and  find  the 
resulting  change  in  pressure  by  equation  (25) .  The  problem  now 
becomes  a  bit  tricky  because  we  can  no  longer  isolate  a  single  flow 
and  solve  for  it  explicitly,  even  if  we  were  to  use  a  quadratic  as 
before.  We  can  get  around  this  problem  by  assuming  that  the  factor 
f  *G  is  very  small.   Then,  we  can  solve  for  Qc  by 


A  (53) 

M  F 

We  can  now  get  Qm  from  equation  (46)  ,  and  solve  the  quadratic  of 
equation  (47)  for  Qr: 

m     -X±JX2-4Y(XQm-AP)  (54) 

2Y 


0r 


Once  Qr  is  known,  Q  can  be  solved  for.  We  check  our 
assumption  of  Qt  by  solving  for  change  in  pressure  in  equation 
(30)  ,  iterating  on  Qt  until  the  two  values  of  change  in  pressure 
are  equal.  Again,  these  calculations  are  done  for  us  on  the 
spreadsheet. 

Finally,  we  check  our  assumption  that  the  quantity  fs*G  is 
negligible  by  solving  for  change  in  pressure  in  equation  (52) ,  and 
comparing  with  our  answer.  By  doing  so,  we  find  that  the 
assumption  was  reasonable  (solution  error  <  2%) .  If,  however,  we 
want  to  "fine-tune"  the  answer,  we  can  now  calculate  the  factor 
fs*G,  plug  it  into  equation  (52) ,  and  solve  the  remaining  equations 
as  before.  As  it  turns  out,  the  factor  fs*G  does  not  vary  much 
with  changing  Qs  over  the  range  of  values  in  which  we  are 
interested;  we  can  therefore  assume  a  constant  value  based  on  a 
single  calculation,  giving  an  improved  (though  still  slightly 
inaccurate)  estimation  of  Qs.  My  calculations  gave  a  value  of  fs*G 
=  .03.   Equation  (53)  then  becomes 


pi  (55) 

\  F+  .03 


Case  #  4 i   Bypass  valve  in  place,  turbulent  flow. 

For  this  problem,  all  equations  of  case  #  3  remain  the  same, 
with  the  exception  of  equation  (47) ,  which  must  now  be  replaced  by 
an  expression  for  turbulent  flow  losses.  The  method  is  the  same  as 
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was  used  in  case  number  two;  in  fact,  the  equations  are  the  same: 

P1   =  WfVJ  +  YQJ  (56) 

Vf  =  [°'*°'\       1       \  (57) 
1      [   40af  j\  60000/ 


We  solve  for  change  in  pressure  and  Qs  using  equations  (25)  &  (55)  , 
and  for  Qm  using  equation  (46)  .  Using  Q  =  Qm-Qr/  we  substitute  into 
equation  (30)  and  get 

P,  =  2Z(Q„-QZ)  -YQr  (58) 

Solving  the  quadratic,  we  obtain 

_  -2Z±JtZ2+*Y(2ZQa-~Zp)  (59) 

r  =  2  Y 

and  solve  for  Qr.  We  can  then  obtain  Q  and  check  our  estimation 
of  Qt  by  solving  for  P1  in  equation  (30) ,  again  iterating  until  the 
values  match. 

Results:  Attempts  to  simulate  the  membrane  flow 
characteristics  obtained  in  the  lab  with  the  hydraulic  model  were 
hampered  because  of  a  lack  of  information  concerning  the  pumps  used 
in  lab  experiments.  Without  knowledge  of  the  pump  curve  for  the 
pump  used,  an  exact  comparison  of  data  could  not  be  made.  This 
obstacle  was  partially  overcome  by  deriving  a  pump  curve  equation 
from  known  data  for  a  pump  which  was  not  used  in  the  experiments. 
The  hydraulic  model  was  then  solved,  and  the  resulting  bypass  flow 
adjusted  until  the  membrane  cartridge  influent  flow  of  the 
hydraulic  model  matched  that  of  the  experiment.  The  bypass  flow 
for  the  model  was  controlled  by  varying  the  K  value  of  the  bypass 
valve.  Once  this  match  was  achieved,  a  comparison  of  such 
parameters  as  trans-membrane  pressure,  ratio  of  permeate  to 
membrane  flow,  and  fibre  flow  velocity  could  be  made. 

Table  (4)  and  table  (5)  show  a  comparison  of  experimental  and 
simulated  membrane  parameters  for  two  different  experiments.  Each 
parameter  is  also  plotted,  as  shown  in  figures  (5-12) .  The  model 
responds  favorably  in  most  cases;  Experiment  number  two  shows 
results  which  vary  considerably  in  cases  for  which  inlet  pressures 
which  differ  from  the  norm.  I  would  speculate  that  the  discrepancy 
lies  in  the  possibility  that  the  bypass  flow  was  altered  during  the 
time  that  these  data  were  collected. 
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TABLE  4 
EVALUATE  BOTH  MODELS  FOR  EXPERIMENT  #1 : 

TEST  CONDITIONS 


TEST  # 
1 

1     c 
|1.57E+11 

1.74 

0.398 

Qr 

1.22 

Qt 
1.618 

Qm 

1.618 

Qp/Qm 

0.246 

Pi 
20 

6 

2 

|1.87E+11 

1.66 

0.387 

1.16 

1.547 

1.547 

0.250 

26 

2 

3 

|1.92E+11 

1.75 

0.381 

1.24 

1.621 

1.621 

0.235 

27 

4 

|1.82E+11 

1.67 

0.402 

1.16 

1.562 

1.562 

0.257 

26 

5 

5 

|  1.9E+11 

1.69 

0.393 

1.18 

1.573 

1.573 

0.250 

27 

1 

AVG: 

|1.82E+11 

1.70 

0.392 

1.19 

1.584 

1.584 

0.247 

25 

48 

LAMINAR  FLOW  MODEL 

TEST  # 
1 

1        c 

|1.57E+11 

Vf 
1.680 

Qp 
0.428 

Qr 

1.155 

Qt 
2.152 

Qm 

1.584 

Qp/Qm 

0.271 

Pi 

24 

07 

2 

|1.87E+11 

1.684 

0.376 

1.185 

2.142 

1.561 

0.241 

25 

11 

3 

|1.92E+11 

1.685 

0.369 

1.189 

2.140 

1.558 

0.237 

25 

26 

4 

|1.82E+11 

1.684 

0.384 

1.180 

2.143 

1.564 

0.246 

24 

97 

5 

|  1.9E+11 

1.684 

0.372 

1.187 

2.141 

1.559 

0.239 

25 

20 

AVG: 

|1.82E+11 

1.684 

0.384 

1.180 

2.143 

1.564 

0.246 

24 

96 

TURBULENT 

FLOW  MODEL 

TEST  # 

1 

1        c 

|1.57E+11 

Vf 
1.490 

Qp 
0.451 

Qr 

0.989 

Qt 
2.082 

Qm 

1.440 

Qp/Qm 

0.313 

Pi 
30 

75 

2 

|1.87E+11 

1.497 

0.396 

1.022 

2.071 

1.418 

0.279 

31 

78 

3 

|1.92E+11 

1.498 

0.388 

1.027 

2.070 

1.415 

0.274 

31 

93 

4 

|1.82E+11 

1.496 

0.404 

1.017 

2.073 

1.421 

0.284 

31 

63 

5 

|  1.9E+11 

1.497 

0.391 

1.025 

2.070 

1.416 

0.276 

31 

87 

AVG: 

I1.82E+11 

1.496 

0.404 

1.017 

2.073 

1.421 

0.284 

31 

63 

TABLE  5 
EVALUATE  BOTH  MODELS  FOR  EXPERIMENT  #  2: 


TEST  CONDITIONS 


"1 

TEST  #: 

C 

vf 

QP 

Qr 

Qt 

Qm 

Qp/Qm 

Pi 

1 

1.89E+11 

0.874 

0.145 

0.640 

0.785 

0.785 

0.185 

9.64 

2 

1.95E+11 

1.296 

0.213 

0.950 

1.163 

1.163 

0.183 

14.36 

3 

2.00E+11 

1.085 

0.179 

0.795 

0.973 

0.973 

0.183 

11.77 

4 

2.04E+11 

0.704 

0.126 

0.511 

0.637 

0.637 

0.198 

7.71 

5 

2.10E+11 

1.088 

0.153 

0.811 

0.963 

0.963 

0.159 

9.80 

6 

2.11E+11 

0.888 

0.149 

0.649 

0.798 

0.798 

0.187 

9.40 

7 

2.14E+11 

0.766 

0.288 

0.480 

0.768 

0.768 

0.375 

19.58 

8 

2.16E+11 

0.956 

0.078 

0.740 

0.818 

0.818 

0.096 

4.64 

AVG: 

2.05E+11 

0.957 

0.166 

0.697 

0.863 

0.863 

0.182 

10.86 
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LAMINAR  FLOW  MODEL 


TEST  #: 

C 

vf 

Qp 

Qr 

Qt 

Q* 

Qp/Qn, 

Pi 

1 

1 

89E+11 

0.940 

0.144 

0.694 

2.286 

0.838 

0.172 

10 

45 

2 

1 

95E+11 

0.939 

0.140 

0.696 

2.285 

0.836 

0.168 

10 

48 

3 

2 

00E+11 

0.939 

0.136 

0.697 

2.285 

0.833 

0.164 

10 

51 

4 

2 

04E+11 

0.938 

0.134 

0.698 

2.285 

0.832 

0.162 

10 

52 

5 

2 

10E+11 

0.938 

0.131 

0.699 

2.285 

0.830 

0.158 

10 

55 

6 

2 

11E+11 

0.938 

0.130 

0.699 

2.285 

0.830 

0.157 

10 

56 

7 

2 

14E+11 

0.937 

0.129 

0.700 

2.284 

0.828 

0.155 

10 

57 

8 

2 

16E+11 

0.937 

0.128 

0.700 

2.284 

0.828 

0.154 

10 

58 

AVG: 

2 

05E+11 

0.938 

0.134 

0.698 

2.285 

0.832 

0.161 

10 

53 

TURBULENT  FLOW  MODEL 


"1 

TEST  # 

C 

vf 

Qp 

Qr 

Qt 

Qm 

Qp/Qn, 

Pi 

1 

1.89E+11 

0.850 

0.142 

0.622 

2.277 

0.764 

0.186 

11.42 

2 

1.95E+11 

0.850 

0.138 

0.624 

2.277 

0.762 

0.181 

11.45 

3 

2.00E+11 

0.850 

0.135 

0.625 

2.276 

0.760 

0.177 

11.47 

4 

2.04E+11 

0.849 

0.133 

0.626 

2.276 

0.759 

0.175 

11.49 

5 

2.10E+11 

0.849 

0.129 

0.628 

2.276 

0.757 

0.171 

11.51 

6 

2.11E+11 

0.850 

0.129 

0.629 

2.276 

0.757 

0.170 

11.51 

7 

2.14E+11 

0.849 

0.127 

0.629 

2.276 

0.756 

0.168 

11.53 

8 

2.16E+11 

0.849 

0.126 

0.629 

2.276 

0.755 

0.167 

11.53 

AVG: 

2.05E+11 

0.850 

0.132 

0.626 

2.276 

0.759 

0.174 

11.49 

Lamir 

iar  vs.    Turbulent  modeling 

:      Both  mode 

Is   seem  to   s 

expe:.  imental  data  reasonably  well.  This  tends  to  support  the  idea 
that  the  membrane  fibers  operate  in  the  transition  region  between 
laminar  and  turbulent  flows.  It  should  be  noted  that,  by 
adjusting  the  K  value  of  the  bypass  valve  in  the  hydraulic  model, 
one  could  cause  one  model  to  perform  better  than  the  other.  This 
is  possible  because  of  the  turbulent  model's  increased  resistance 
to  fibre  flow.  This  increased  resistance  tends  to  increase  the 
cartridge  inlet  pressure,  which  can  be  accurately  simulated  by 
increasing  the  K  value  of  the  bypass  valve  (this  in  effect  "closes" 
the  bypass  valve,  thereby  increasing  inlet  pressure  as  well  as 
increasing  flow  to  the  membrane) . 

Discussion:  The  hydraulic  model  as  written  is  an  effective 
simulation  of  hollow  fibre  membrane  flow.  The  model  also  contains 
the  flexibility  to  simulate  fouling  of  the  membrane,  as  well  as 
hollow  fibre  membranes  of  varying  composition  (see  appendix  A  for 
a  description  of  how  the  spreadsheet  operates).  The  model's 
limitation  is  that  it  is  static  in  nature;  although  a  dynamic 
system  cannot  be  simulated  directly,  it  can  be  configured  to 
simulate  incremental  changes  in  time.  This  is  done  by  making  the 
membrane  constant  a  function  of  time,  solving  for  the  membrane 
constant  for  different  values  of  time,  then  imputing  that  value 
into  the  model  and  solving  the  resulting  iteration.  The  membrane 
constant  may  be  viewed  in  this  regard  as  the  degree  to  which  the 
membrane  is  fouled.  If,  for  example,  one  believes  the  membrane  to 
foul  exponentially,   one  could  make  the  membrane  constant  an 
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exponential  function  of  time,  then  solve  the  model  for  various 
values  of  time.  This,  as  we  shall  see,  is  the  approach  used  in  the 
section  on  economic  optimization. 

Further  research:  Several  assumptions  had  to  be  made  in  the 
development  of  these  results;  specifically,  the  use  of  a  pump  curve 
which  was  not  representative  of  the  pumps  used  in  the  actual 
experiments,  and  a  speculation  of  the  values  of  the  bypass  valve 
constant  K  and  bypass  flow  rate.  Further  research  experiments 
should  be  designed  to  specifically  measure  bypass  flow  rate  and 
change  in  pressure  across  the  bypass  valve,  as  well  as  permeate  and 
retentate  flows  and  trans-membrane  pressure.  Further  research 
experiments  should  also  include  pumps  for  which  the  characteristic 
pump  curve  is  known.  This  information  is  essential  in  improving 
the  accuracy  of  the  hydraulic  model,  as  well  as  in  providing 
further  insight  to  the  preference  of  the  laminar  flow  model  over 
the  turbulent  flow  model.  Finally,  as  pointed  out  in  Clark,  1991, 
an  effort  should  be  made  to  incorporate  an  irreversible  fouling 
term  into  the  model.  This  could  most  easily  be  done  by  making  the 
membrane  constant  an  increasing  function  of  permeate  flow.  The 
function  which  the  term  would  take  would  require  considerable 
research  into  the  fouling  mechanisms  of  the  membrane  being  studied. 
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Figure    (5) 
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Figure    (6) 
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Figure    (7) 
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Figure    (3) 
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Figure    (9) 
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Figure    (10) 
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Figure    (11) 
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Figure    (12) 


ECONOMIC  OPTIMIZATION 


The  cost  of  UF  operations  has  been  discussed  in  several  EPA 
reports  and  in  other  research  (Gumerman  et  al,  1985,  Taylor  et  al, 
1989,  and  Laine  et  al,  1991).  Gumerman  et  al  (1985)  developed 
capital  cost  estimates  in  consideration  of  treatment  plant 
construction  planning.  Taylor  et  al  (1989)  developed  operating 
cost  estimates,  to  include  membrane  replacement,  for  life  cycle 
costs  associated  with  the  bench-scale  use  of  spiral  wound 
membranes.  In  their  analysis,  Taylor  et  al  used  the  mass  transfer 
theory  model  to  estimate  permeate  flux.  Laine  et  al  (1991) 
related  operating  costs  to  UF  pump  energy  costs,  without  regard  for 
membrane  replacement.  This  report  will  attempt  to  optimize 
operating  costs  of  a  UF  operation.  Amortization  of  capital  costs 
will  not  be  considered  in  this  analysis. 

Because  of  fouling  of  the  UF  membrane  during  the  normal  course 
of  operations,  it  is  necessary  to  backflush  the  membrane 
periodically.  This  backflush  serves  both  to  increase  the  permeate 
flux  and  to  reduce  the  pressure  head  that  the  pump  must  work 
against,  thus  making  the  process  less  expensive.  Back  flushing, 
however,  requires  a  source  of  clean  water  and,  of  course,  energy. 
The  clean  water  source  is  simply  the  permeate  previously  generated. 
We  will  show  that  there  is  an  optimal  cycle  time  (one  cycle 
includes  one  UF  period  and  one  backflush  period)  at  which  the 
system  should  operate;  that  is  to  say,  a  unique  cycle  time  exists 
for  which  the  volume  of  permeate  produced  per  kilowatt-hour  of 
energy  expended  is  maximized. 

Our  intention  is  to  quantify  both  the  permeate  produced  and 
energy  consumed  as  functions  of  time.  We  may  then  perform  some 
basic  calculus  to  find  the  maximum.  The  energy  required  of  any 
pump  operating  at  a  given  flow  rate  and  pressure  head  is  given  by 

Power  =  pgQhp  (60) 

Assuming  a  pump  efficiency  of  80%,  the  energy  required  to  produce 

flow  Q  at  head  h  is 

p 

Power  -    pgQh*  (61) 

.8 


Substituting  the  relationship  AP=pgh  ,  we  get  the  following: 

Power   =  1.25Q&P  (62) 

Applying  conversion  factors  for  PSI  to  N/M2  and  M3/sec  to  1/min,  the 
end  result  is 

Power   =  .1447  CAP  (63) 

where  power  is  measured  in  watts.  We  therefore  may  find  the  power 
delivered  to  either  the  UF  pump  or  the  backflush  pump  if  we  know 
both  the  discharge  pressure  and  the  flow  delivered. 

Now  that  the  power  relationship  has  been  defined,  we  need  to 

28 


express  both  power  and  permeate  flow  a  functions  of  time.  Both 
power  and  permeate  flow  have  two  components;  one  for  UF,  and  one 
for  backflush.  At  this  point,  we  need  to  make  some  rather  bold 
assumptions  about  the  system  hydraulics.  Our  first  assumption  is 
that  the  UF  pump  is  operating  in  the  high  flow,  low  pressure  part 
of  the  pump  curve,  as  indicated  in  figure  (13) . 


AP 
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Q  (1/min) 


r-rr 


Figure  (13) 


This  conforms  with  the  current  laboratory  set  up,  and  allows 
us  to  make  the  assumption  that  Qt,  the  total  UF  discharge,  is 
independent  of  pressure  (over  the  range  of  values  in  which  we  are 
interested) .  Although  this  is  an  approximation,  it  is  reasonable; 
the  hydraulic  model  shows  that  Qt  changes  only  about  2%  over  the 
range  of  pressures  in  which  we  are  interested.  Our  last,  and 
perhaps  most  inaccurate  assumption,  is  defining  how  Q  varies  with 
time.  In  reality,  Q  is  a  function  of  the  rate  at  which  the 
membrane  is  fouled.  Lacking  sufficient  experimental  data,  we  will 
use  an  exponential  decay  function  accepted  by  many  researchers 
(Cheryan,  1986,  p. 174): 


Op  '  oPot-™ 


(64) 


Experimental  data  shows  a  linear  decay  of  the  permeate  flux  for  the 
first  15  minutes  of  UF  operation  (Heneghan,  1991,  Adham,  1991)  , 
dropping  off  to  about  75%  of  the  original  permeate  flux.  The 
literature  indicates  that  the  decay  rate  slows  considerably  after 
this  initial  decline  (Cheryan,  1986,  p. 172;  Fell,  1990).  The 
coefficient  selected  in  equation  (64)  approximates  these  results, 
and  plots  as  shown  in  figure  (14)  ;  it  is  characterized  by  a 
relatively  linear  decay  for  the  first  15  minutes  of  the  cycle,  then 
exponentially  decaying  to  approximately  75%  of  the  initial  value. 
The  curve  generated  by  equation  (64)  is  a  "best  guess,"  and  should 
be  improved  upon  by  comparison  with  more  extensive  experimental 
data.  Integrating  equation  (64) ,  we  get  permeate  volume  as  a 
function  of  time: 


Vp(t)  =  1.0810  c"5 


(65) 


Now  having  an  expression  for  permeate  flow  as  a  function  of  time, 
we  can  determine  permeate  flow  at  any  point  in  time  by  simple 
calculation.  We  can  then  take  these  generated  values,  plug  them 
back  into  the  hydraulic  model,  and  solve  for  the  membrane  constant, 
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Figure    (14) 


C,  which  corresponds  to  the  desired  flow.  Each  flow  will  also 
generate  a  unique  discharge  pressure  (AP) ,  which  we  can  plot 
against  the  corresponding  times  for  which  the  permeate  flows  were 
taken.  Finally,  we  can  fit  a  mathematical  expression  to  this 
curve.  For  the  data  in  this  example,  the  expression  for  AP  as  a 
function  of  time  is: 

p  =  P0*5.66  (l-f-1)  (66> 

Figure  (15)  shows  the  resulting  curve.  Note  that  this  curve  is 
based  on  the  data  generated  by  the  spreadsheet  from  equation  (64), 
and  should  be  improved  upon  by  comparison  with  more  extensive 
experimental  data. 

Now,  substituting  equation  (66)  into  the  power  equation  yields 

Power   =  .  1447  0(Po  +  5.66  (l-f1)  )  (67) 

Integrating  this  expression  over  time  results  in  an  equation  for 
the  energy  delivered  to  the  pump  (note  unit  conversions  from  watts 
to  kilowatts  and  from  minutes  to  hours) : 

kW-H(C)    =  2.41xlO-6Qc  [t(P0+5.66)  -6  .29  f9]  (68) 

Our  goal  is  to  optimize  the  ratio  of  Permeate  volume  to 
kilowatt-hours  consumed.  In  this  regard,  we  must  now  account  for 
the  power  required  to  drive  the  backflush  pump  as  well  as  the 
permeate  volume  required  for  the  backflush.  We  do  so  as  follows: 
We  will  make  the  assumption  that  the  membrane  module  experiences  no 
irreversible  fouling.  We  will  further  assume  that  a  one  minute 
backflush  will  achieve  a  100%  removal  of  particulate  matter  from 
the  membrane  (or,  100%  recovery  of  the  initial  permeate  flux) . 
Experimental  data  for  the  laboratory  set  up  we  have  used  indicated 
a  backflush  pump  discharge  pressure  of  about  33  PSI,  with  a  flow 
rate  of  .35  1/min.  Assuming  a  constant  backflush  discharge 
pressure,  we  obtain  the  following  results: 

kWH  =  l.AlxlO^QLPLt  =  2.41x10-*  (.35)  (33)  (1) 

kW-H  =  2.11xl0-bkW-H  (69) 

Vol.    =     ,35J  x lmin  =   .35  2 
min 

It  is  important  to  note  that  we  are  assuming  that  a  one  minute 
backflush  will  suffice  regardless  of  the  length  of  the  period  of  UF 
prior  to  the  backflush.  Intuitively,  we  would  suspect  that  this  is 
not  the  case,  and  that  the  length  of  the  backflush  would  in 
actuality  be  a  function  of  the  length  of  the  period  of  UF.  For  the 
moment,  we  will  proceed  as  planned,  and  see  if  our  assumption  is 
accurate. 

From  the  above  equations  we  can  now  predict  the  permeate 
volume  generated  and  the  energy  consumed  in  one  cycle  (one  cycle  = 
one  UF  and  one  backflush  period) .  Combining  equations  (65) ,  (68) , 
and  (69) ,  we  get  the  desired  equations  for  energy  and  permeate 
volume: 
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Figure    (15) 


kWH{t)    =  2.41XJ0-SCC  [t(P0+5.66)  -6.29  f»]  +2.77XJ0-5  (70) 

Vp(t)    =  1.081CPot-925-  .35  (71) 

We  now  must  determine  the  time  for  which  the  ratio  of  power  to 
volume  is  at  a  maximum.  This  can  be  most  simply  done  by  dividing 
equation  (71)  by  equation  (70)  ,  then  take  the  first  derivative  with 
respect  to  time.  Setting  the  resulting  expression  to  zero  and 
solving  for  time  will  provide  us  with  the  maximum.  The  form  of  the 
derivative  is 


d 
dt 


Vp(t)    1  _    1    (dVp(t))  VO   ldkW-H(t)\  (72) 

{kW-H(t))~     kWH(t){       dt       )       (kW-H(t))2\  dt         ) 


Once  the  derivative  is  taken,  we  set  the  resulting  expression  equal 
to  zero  and  solve  for  time.  The  mathematics  of  doing  this  are 
somewhat  complicated,  so  instead  we  will  use  the  spreadsheet  and 
solve  by  iteration.  The  results  are  shown  in  table  (6)  .  A 
graphical  solution  is  also  possible.  By  dividing  equation  (71) 
into  equation  (70)  and  solving  the  resulting  expression  for  various 
values  of  time  (table  7) ,  we  can  plot  the  results  as  indicated  in 
figure  (16) ,  which,  as  we  expect,  provides  the  same  conclusion  as 
the  calculus  solution. 


TABLE  6 
CALCULUS  SOLUTION  OF  MAXIMIZATION  OF  t  FOR  EXPONENTIAL  DECAY: 


Qt 

AP 


t  - 


0.391 

VP  ~ 

4.185 

2.119 

dVp/dt  - 

0.323 

28.181 

kWH  - 

1.95E-03 

dkWH/dt  - 

0.00015 

13.007 

d/dt(VpAWH)- 
TABLE  7 

0.00176 

GRAPHICAL  SOLUTION 
OPTIMIZATION  OF  PERMEATE  FLOW  PER  KILOWATT  HOUR 


t 

vP 

kW-H 

Vp  per  kW-H 

(min) 

(1/min) 

4 

1.17 

6.0708E-04 

1933.4 

6 

1.87 

9.0346E-04 

2066.7 

8 

2.54 

1.2015E-03 

2116.6 

10 

3.21 

1.5007E-03 

2136.6 

12 

3.86 

1.8008E-03 

2143.2 

14 

4.50 

2.1017E-03 

2143.4  *** 

16 

5.14 

2.4033E-03 

2140.0 

18 

5.78 

2.7053E-03 

2134.8 

20 

6.40 

3.0079E-03 

2128.5 

22 

7.02 

3.3109E-03 

2121.7 
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Figure    (16) 


24 

7.64 

3.6143E-03 

2114.6 

26 

8.26 

3.9180E-03 

2107.5 

28 

8.87 

4.2221E-03 

2100.3 

30 

9.48 

4.5264E-03 

2093.3 

32 

10.08 

4.8310E-03 

2086.4 

34 

10.68 

5.1359E-03 

2079.7 

36 

11.28 

5.4410E-03 

2073.1 

38 

11.88 

5.7464E-03 

2066.8 

40 

12.47 

6.0519E-03 

2060.6 

Discussion/further  research:  The  procedure  outlined  above  is  an 
effective  procedure  to  optimize  UF  operations.  Its  weakness  lies 
in  the  assumptions  that  were  required,  specifically,  the  response 
of  the  UF  system  to  long-term  fouling,  and  the  resulting  effect  on 
inlet  pressure.  Further  research  should  include  an  experiment 
which  measures  permeate  flow  rate  and  membrane  cartridge  inlet 
pressure  over  long  periods  of  time  in  which  there  is  no  back  flush 
operations.  Attention  should  also  be  given  to  the  amount  of  time 
required  for  complete  cleaning  of  the  membrane  as  a  function  of 
cycle  time  (recall,  we  assumed  a  one  minute  back  flush  would 
suffice  regardless  of  the  length  of  the  UF  portion  of  the  cycle) . 
An  evaluation  of  the  actual  time  required  for  backflush  as  the  UF 
period  becomes  longer  would  be  beneficial  in  further  improving  the 
accuracy  of  the  optimization.  Finally,  the  optimization  procedure 
does  not  account  for  the  size  of  the  UF  pump.  The  pump  we  used 
operated  in  the  high  flow,  low  pressure  region  of  the  pump 
operating  curve.  This  suggests  that  the  pump  was  too  large  for  the 
system  requirements,  resulting  in  considerable  bypass  flow;  as  a 
result,  one  could  make  a  case  for  a  smaller  pump  which  operates  at 
a  lower  cost  per  kilowatt  hour,  and  which  has  a  lower  capital  cost 
(an  issue  not  previously  discussed,  but  outside  the  scope  of  this 
report) .  The  trade  off  in  this  case  is  that,  although  operating 
costs  are  initially  lower,  Qt  would  no  longer  be  independent  of 
membrane  cartridge  inlet  pressure.  This  effect  would  necessitate 
a  change  in  the  equation  for  pump  power  (equation  (67)),  which  must 
then  account  for  a  change  in  Q  as  a  function  of  time.  The 
corrections  could  be  made  with  relative  ease,  and  the  rest  of  the 
procedure  followed  as  indicated. 
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APPENDIX  1 
SPREADSHEET  OPERATION 


The  spreadsheets  developed  in  conjunction  with  this  report  were  formulated 
with  QUATTRO  PRO .  version  3.0.  In  all,  three  different  spreadsheets  were 
written;  one  for  the  model  in  which  the  bypass  valve  was  not  utilized,  one  for 
which  the  bypass  line  was  utilized,  and  one  for  the  economic  optimization.  The 
first  two  were  designed  specifically  for  future  manipulation,  including 
evaluation  of  different  pumps,  evaluation  of  the  effect  of  fouling,  and  for  ease 
of  improvement  as  more  experimental  data  becomes  available.  The  third  spread 
sheet  is  little  more  than  a  scratch  sheet  used  to  solve  the  economic 
optimization;  because  of  its  relatively  simple  format,  it  will  not  be  further 
discussed.  The  intention  of  the  rest  of  this  appendix  is  to  familiarize  the 
reader  with  the  use  and  application  of  the  first  two  spreadsheets.  Details  on 
the  use  of  QUATTRO  PRO  may  be  found  in  the  user's  guide. 

The  spreadsheet  for  the  hydraulic  model  without  the  bypass  line  is  labeled 
"M0DEL#3 . 1 . "  The  spreadsheet  for  the  model  with  the  bypass  line  is  labeled 
"M0DEL#4 . 1 . "  The  spreadsheets  are  set  up  identically,  and  are  operated  in  the 
same  manner.  For  explanation  purposes,  we  shall  use  M0DEL#4 . 1 .  Figure  A-l 
shows  the  layout  of  the  operating  portion  of  the  spreadsheet.  The  spreadsheet 
is  designed  to  solve  for  a  number  of  system  parameters  based  on  an  input  value 
of  the  membrane  constant,  C.  Once  a  value  of  C  is  entered  into  the  spreadsheet, 
a  non-zero  value  of  "Check  PI"  will  appear  in  the  input.  Qt  is  then  iterated 
upon  until  "Check  PI"  is  equal  to  zero.  This  iteration  is  most  efficiently 
accomplished  by  using  the  /Tools/Solve  For  command,  which  can  run  up  to  99 
iterations.  The  iteration  must  be  performed  for  both  the  laminar  and  the 
turbulent  models.  Once  the  iteration  is  complete,  the  LAMINAR/TURBULENT  FLOW 
MODEL  blocks  will  indicate  the  output  values  of  the  parameters  being  evaluated. 

In  the  event  that  a  series  of  output  values  are  to  be  evaluated  (or 
graphed) ,  such  data  may  be  recorded  in  the  following  manner  (refer  to  figure  A- 
2):  Use  the  [ctrl]  C  command  to  copy  the  lower  half  of  the  output  template  onto 
the  upper  half  as  shown  in  figure  A-2 .  This  will  replace  the  old  data  in  the 
upper  half  of  the  template  with  the  cell  addresses  of  the  output  parameters  in 
figure  A-l.  As  each  set  of  output  values  is  generated,  use  the  /Edit/Values 
command  to  replace  the  cell  addresses  in  the  upper  half  of  the  template  with  the 
actual  numbers  shown.  This  procedure  in  effect  preserves  the  cell  address  matrix 
for  future  use.  Note  that  the  values  of  the  membrane  constant  indicated  on  the 
template  are  for  reference  purposes  only.  Because  they  are  not  required  for  any 
calculations,  they  may  be  replaced  with  the  values  used  to  generate  the  new 
output.  The  reader  is  referred  to  the  QUATTRO  PRO  user's  guide  for  information 
concerning  the  generation  of  graphs. 

The  input  window  in  figure  A-l  also  has  a  cell  for  input  of  the  bypass 
valve  constant  K.  This  value  may  be  increased  ("closing")  or  decreased 
("opening")  as  required  by  the  operator.  Once  experimental  data  is  available  to 
effectively  estimate  K,  it  should  no  longer  require  adjustment. 

The  spreadsheet  is  also  designed  to  allow  for  changes  in  the  lab 
configuration  or  the  membrane  module  used.  Figure  A-3  shows  the  input  window  for 
the  parameters  of  interest.  Such  parameters  as  tubing  length  and  diameter,  pump 
curve  constants,  viscosity,  density,  and  membrane  parameters  may  be  changed. 
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EVALUATE  BOTH  MODELS  FDR  EXPERIMENT  #1: 


TEST  CONDITDNS 


TEST 
1 
2 
3 
4 
5 
AVG: 


#: 


c 

Vf 

Op 

Or 

Qt 

Qm 

Qp/Qm 

P1 

1.569E+11 

1.74 

0.398 

1.22 

1.616 

1.618 

0.24B 

20.6 

1.874E+11 

1.BB 

0.387 

1.16 

1.547 

1.547 

D250 

2B2 

1.91BE+11 

1.75 

0.3B1 

1.24 

1.B21 

1.B21 

0235 

27 

1.B21E+11 

1.B7 

0.402 

1.16 

1.5B2 

1.5B2 

0.257 

26.5 

1.9Q5E+11 

1.69 

0.393 

1.16 

1.573 

1.573 

0.250 

27.1 

1.B1BE+11 

1.7D 

0.392 

1.192 

1.584 

1.584 

0.247 

25.48 

TEST 
1 
2 
3 
4 
5 
AVG: 


ft 


LAMNAR  FLOW  MODEL 

c 

Vf 

Op 

Or 

Dt 

Qm 

Qp/Qm 

P1 

1.569E+11 

1.BB0 

0.42B 

1.155 

2.152 

1.584 

0.271 

24.073 

1.B74E+11 

1.BB4 

0.376 

1.185 

2.142 

1.5B1 

0.241 

25114 

1.91BE+11 

1.BB5 

0.369 

1.189 

2.140 

1.55B 

0.237 

25.264 

1.B21E+11 

1.BB4 

0.3B4 

1.180 

2.143 

1.564 

0.246 

24.956 

1.9D5E+11 

1.BB4 

0.372 

1.167 

2.141 

1.559 

0.239 

25.205 

1.B1BE+11 

1.BB4 

0.384 

1.180 

2.143 

1.564 

0.246 

24.956 

TEST 
1 
2 
3 
4 
5 
AVG: 


ft 


TURBULENT  FLOW  MODEL 

c 

Vf 

Qp 

Or 

Qt 

Qm 

Qp/Qm 

P1 

1.569E+11 

1.490 

0.451 

0.989 

2.082 

1.440 

0.313 

30.746 

1.674E+11 

1.497 

0.396 

1.022 

2.071 

1.418 

0.279 

31.7B3 

1.91BE+11 

1.49B 

0.388 

1.027 

2.070 

1.415 

0.274 

31.929 

1.B21E+11 

1.49B 

0.404 

1.017 

2.073 

1.421 

0.264 

31.627 

1.905E+11 

1.497 

0.391 

1.025 

2.070 

1.416 

0.276 

31.874 

1.81BE+11 

1.496 

0.404 

1.017 

2.073 

1.421 

0.284 

31.627 

STORE  TEST  RESULTS  MATRDK 


LAMNAR  FLOW  MODEL 

EST  #: 

C 

Vf 

Qp 

Or 

Qt 

Qm 

Qp/Qm 

P1 

1 

1.569E+11 

+£Q£29 

+£M$26 

+$N£26 

+^$12 

+£Q£26 

+£R£26 

+$N£34 

2 

1.B74E+11 

+£Q£29 

+£M£26 

+SNS26 

+£R$12 

+£Q£26 

+ffi£26 

f$N$34 

3 

1.91BE+11 

+£Q$29 

+£M£26 

+£N$26 

+$R£12 

+$Q£26 

+£R£26 

+SNK34 

4 

1.B21E+11 

+£Q$29 

+$M£26 

+£N$26 

+$R£12 

+£Q£26 

+£R$26 

+SNK34 

5 

1.905E+11 

+$Q£29 

+£M$26 

+$N$26 

+$R$12 

+$Q$26 

+£R£26 

+SMS34 

AVG: 

1.B1BE+11 

+£Q$29 

+$M£26 

+£N$26 

+£R£12 

+£Q£26 

+$R£26 

+$N£34 

TURBULENT  FLOW  MODEL 

EST  #: 

C 

Vf 

Qp 

Or 

Qt 

Qm 

Qp/Qm 

dP 

1 

1.5B9E+11 

+$W£3 

+$S$26 

+$T£26 

+$R$1B 

+$W£2 

+SXS26 

+$T£34 

2 

1.B74E+11 

+£W£3 

+$S$26 

+$T£26 

+SRS1B 

+$W£2 

+£X$26 

+£T£34 

3 

1.91BE+11 

+$W£3 

+$S$26 

+£T£26 

+SRS1B 

+£W$2 

+£X£26 

+£T£34 

4 

1.B21E+11 

+SWS3 

+SSS26 

+$T£26 

+£R£1B 

+£W£2 

+£X£26 

+£T£34 

5 

1.905E+11 

+$W£3 

+£S£26 

+£T$26 

+$R$18 

+$W$2 

+£X£26 

+£T£34 

AVG: 

1.81BE+11 

+£W£3 

+£S$26 

+£T$26 

38 

+£R£18 

+£W$2 

+5XS26 

+£T£34 

Figure    (A-2) 

SYSTEM  PARAMETERS: 

Q's  ARE  N  l/mh    dP  B  ISI  psig,    AREAS  ARE  N  rrf2 

LENGTHS  AND  DIAMETERS  ARE  N  m 


PARAMETERS: 

BB95 
BDDDD 

PUMP  CURVE:    A  = 

A  = 

-36.943 

B  = 

B  = 

B1.BB335 

C  - 

C  = 

6340662 

MEMBRANE  CONSTANT  = 

C  = 

2.05E+11 

VBCOSfTY  = 

u  = 

9.41  BE-  04 

CARTRIDGE  SURFACE  AREA  = 

Acart  ■- 

--  006 

FBRE  LENGTH  = 

Lf  = 

1.02 

FBRE  DIAMETER  = 

df  = 

9.30E-04 

FBRE  CROSS  SECTDN  = 

af  = 

6.793E-07 

FBRE  ROUGHNESS: 

ef  = 

7.64E-05 

DENSITY  OF  NFLUENT  = 

P  = 

997.4 

BACK  PRESSURE  VALVE  CONST.  . 

:          K     = 

2530 

CROSS  SECTDN  OF  RETENTATE  1 

'    Abp  = 

6.362E-05 

CROSS  SECTDN  OF  NF.  TUBE: 

At  = 

6362E-05 

LENGTH  OF  NF.  TUBNG 

Lt  = 

0.6 

NF.  TUBE  DIAMETER: 

dt  = 

0.009 

dZ  FROM  PUMP  TO  MEMBRANE: 

dZ  = 

0.6 

NF.  TUBE  ROUBHNESS 

e  = 

1.5E-06 

BYPASS  VALVE  K: 

K  = 

1000 

BYPASS  TUBE  LENGTH: 

Ls  = 

1 

CROSS  SECTDN  OF  BYPASS  TUBNG  * 

6.362E-05 

DIAMETER  OF  BYPASS  TUBE: 

ds  = 

0009 

39 


Figure    (A-3) 
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